Optimal transport and von Neumann entropy in an Heisenberg XXZ chain out of equilibrium 
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In this paper we investigate the spin currents and the von Neumann entropy (VNE) of an Heisenberg XXZ 
chain in contact with twisted XY-boundary magnetic reservoirs by means of the Lindblad master equation. 
Exact solutions for the stationary reduced density matrix are explicitly constructed for chains of small sizes by 
using a quantum symmetry operation of the system. These solutions are then used to investigate the optimal 
transport in the chain in terms of the VNE. As a result we show that the maximal spin current always occurs 
in the proximity of extrema of the VNE and for particular choices of parameters (coupling with reservoirs and 
anisotropy) it can exactly coincide with them. In the limit of strong coupling we show that minima of the VNE 
tend to zero, meaning that the maximal transport is achieved in this case with states that are very close to pure 
states. 
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I. INTRODUCTION 

There is presently an increasing amount of interest for the 
understanding of the transport properties in quantum open 
(out of equilibrium) systems [Ql 0] due to possible applica- 
tion in the field of condensed matter physics and of quantum 
computing. In particular, transport properties in anisotropic 
XXZ Heisenberg spin chains in contact with magnetic reser- 
voirs have been extensively investigated in the past years (for 
extended reviews see Ref. IH @1) with a range of alternate 
methods, including exact diagonalizations Jsllg], Bethe-ansatz 
ill-I^, Lagrange multipliers ifTol [Till . Lanczos method lfl2ll . 
quantum Monte Carlo ifisll . etc. On the other hand, it is well 
known that these systems can be also treated by replacing ob- 
servables associated with the reservoirs (or more in general 
with the environment) by means of quantum noise and study- 
ing the remaining degrees of freedoms within the formalism of 
the quantum open systems Within the Markovian approx- 
imation this amounts to model reservoirs by means of Lind- 
blad operators acting at the system ends and to replace the 
Liouville equation for the time evolution of the density ma- 
trix with a master equation of Lindblad type for the reduced 
density matrix acting of the Hilbert space associated to the 
degrees of freedoms of physical interest filial- 

The aim of the present paper is to use the formalism of the 
Lindblad master equation to provide a first characterization of 
the optimal transport in an open spin chain out of equilibrium 
in terms of the von Neumann entropy (VNE). In this respect 
we consider a XXZ Heisenberg chain with twisted XY bound- 
ary reservoirs at the ends which possesses inversion symmetry 
with respect to the middle point of the chain. By means of this 
symmetry we derive a set of algebraic equations which allow 
to compute in exact manner the stationary density matrix ele- 
ments for chains of small sizes. As a result we show that the 
maximal spin current always occurs in proximity of a mini- 
mum of the VNE and for particular choices of parameters it 
can exactly coincide with it. We also show that in the strong 
coupling limit the minima of the VNE approach zero (exactly 
vanishing for infinite coupling), this meaning that the optimal 



transport is achieved in this case with states that are very close 
to pure states. 

We remark that the open XXZ chain with an effective con- 
stant pumping at the chain ends which induce a gradient 
along the z axis |ll4il or with easy-plane magnetizations at 
the boundaries iflSll . have also been investigated. To the best 
of our knowledge, however, no specific link between optimal 
transport and extrema of VNE has been previously reported. 

The paper is organized as follows. In Section II we present 
the model equation and discuss the quantum symmetry of the 
system. We show (in this section and in the Appendix) how 
the symmetry operator can be used to derive a set of alge- 
braic equations which allow to obtain elements of the station- 
ary reduced density matrix in exact manner. In Section III this 
approach is used to investigate the relationship between spin 
current and optimal transport in terms of the VNE for chains 
of small size. In the last Section the main results of the paper 
are briefly summarized. 

II. MODEL EQUATION AND SYMMETRY PROPERTIES 

We consider the quantum Master equation for the open 
XXZ model in the Lindblad form: 

m m 

where p is the reduced density matrix, H is the XXZ Hamil- 
tonian 

N-l 

H=J2 + + ^^l^Ui} (2) 

k=l 

(we have fixed h — \ and used f to denote adjoint operator) 
and Ljn denote Lindblad operators introducing dissipation in 
the system at the ends of the chain. In the following we con- 
sider the case of only two Lindblad operators of the form 

Li = - ia^{\), U = Gy{N) + ia,(iV), (3) 
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with (Ja{i), * = 1, ...N, a = x,y, z, denoting usual Pauli ma- 
trices acting on site i, N being the length of the chain. This 
choice of Lindblad operators models magnetic reservoirs at 
the left [i ~ 1) and right {i = N) end of the chain, with 
fixed polarization along the axes x and y, respectively. Note 
that the operators Li, L2, introduce boundary gradient in the 
polarization along the negative x and positive y— direction 
which allow to sustain a nonzero spin current in the trans- 
verse z-direction which depends in a non trivial manner on 
the anisotropy parameter A (see also |[l5ll . lfl6ll ). Moreover, it 
is easy to check that with this choice of Lindblad operators 
the master equation becomes invariant under the action of the 
symmetry operator 



T = R- P 



(4) 



where R performs in the common cr| eigenbasis an inversion 
of the spin x-direction followed by a rotation in the x — y 
plane, e.g. 



N 



R=Y[ 



k=l 



i 

1 



(5) 



and the operator P reverses the order of the sites j ^ N — 
j + 1, e.g. P{Ai (g) (g) ... (E) An) = An(S) An-i <S) ■■■ «> Ai, 
with Ai, 2x2 matrices. Notice that P is also performing a 
reflection of the chain with respect to its middle point and as 
such, its action is different for N odd or N even, since for N 
odd the central site is unaffected by P while for N even all 
lattice sites are affected (the reflection occurring with respect 
to an inter-site point). Also, one can readily see that T is a 
unitary operator = T^^ satisfying = 1 and 



TL2 = ^iLi, TLi = iL2 



(6) 



Moreover, T commutes with H (it is a symmetry of the closed 
XXZ chain) as well as with the products L]^Lm, m = 1, 2, 

[T, H] = 0, [T, Ll^Lra] = [T, LraLl] = 0, m = 1, 2, (7) 

this making T, for the chosen Lindblad operators ([3]l, a sym- 
metry of the open XXZ chain. 

In terms of stationary density matrix ps, the above relations, 
together with E g. ([p and the assumption of uniqueness for p 
(verifiable as in iflTIl ). imply that [ps, T] = 0, or equivalently: 



TpsT- 



Ps- 



(8) 



This property can be used to obtain exact unique solutions of 
the stationary Lindblad master equation for arbitrary values of 
the parameters A, 7. In this respect, it is worth to note that the 
stationarity of dU together with the fulfillment of Eq. (O and 
the condition Tr{p) = 1, provide a system of algebraic 
equations for matrix elements pij , with d m given by 



IN 



(1 



22Ar-l 

2^) 



for iVodd, 
for N even. 



(9) 



Thus, for example, for the case = 3, one can easily derive 
from Eq. (|8]l that the most general form for the matrix p, com- 
patible with the symmetry T, is the matrix with real diagonal 
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Figure 1: (Color online) Top panels. Stationary spin current J (left) 
and von Neumann entropy (right) as a function of parameters 7, A 
for the open T-Lxxz chain with Lindblad operators ^ and length 
= 3. Bottom panels. Sections of the spin current J (left) and von 
Neumann entropy (right) surfaces in the top panel taken at different 
values of parameter 7. 



elements given by 

Pii = 
P33 = 



P88 
P66 



an, P22 

^33, P55 



and with off-diagonal elements pij 
isfying the relations: 



- Pii 

- P77 

P*,A 



0.22, 
0-551 



a 



(10) 



ibij sat- 



P28 


-Pli, P3i 


iPw, P3S 


-Pie 


Pi5 


= -iP27, Pi6 


= i/023' ^47 


= P25, 


PiS 


= iPl2i P56 = 


= «P377 P5S = 


= -Pl7> 


P67 


= P35 7 Pes ~ 


«Pl3 7 P7S = 


iPi5, 


Pis 


= (1 + z)ai8, 


P2i = (1 + 


i)a2i, 


Pse 


= (1 + i)a36, 


P57 = (1 + 


i)a57, 



(11) 



where the star denotes the complex conjugation and Uij, bij, 
are real numbers (see details in the Appendix). By substitut- 
ing this form of p into Eq. ([T]) and requiring stationarity, one 
gets a system of 32 independent equations for aij , bij which 
admits a unique solution when the normalization condition 
Tr[p] = 2(aii + 022 + 033 + 055) — 1 is imposed. We remark 
that while for interacting closed many body systems a com- 
plete characterization of the reduced density matrix in terms 
of the system symmetries is possible (mainly for systems with 
the permutational invariance ifisll ). very little is known in this 
respect for open quantum systems. 

Although our approach is completely general, the rapid ex- 
ponential growth of the number of algebraic equations to solve 
as N is increased see Eq. quickly restricts the application 
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Figure 2: (Color online) Average spin current vs von Neumann en- 
tropy for different values of the coupling parameter 7. The curves 
are obtained as parametric plots for fixed 7 and A varied in the range 
] — cxD, oo[. The VNE attains it maximum value (same for all curves) 
at J = and for A — ±00. 



of the symmetry method to chains of small sizes. On the other 
hand, qubits systems with few sites can be already of potential 
interest for applications, and in view to the scarcity of exact 
results for open systems, we shall take the opportunity to in- 
vestigate transport properties of these systems in exact manner 
for arbitrary parameters. In this respect we remark that this 
problem has been also investigated in previous papers using 
quantum trajectory method |[l4ll o r direct time integrations of 
the Lindblad master equation flSt l. Our approach, being exact, 
is free of convergence problems of time integrations present 
in these alternate approaches (typically the convergence be- 
comes very slow in the critical region — 1 < A < 1 and in the 
weak and strong coupling limits). 

In addition, the steady state density matrices for positive 
and negative values of the anysotropy parameter A are con- 
nected via relations [U6|] 

PsiN, -A) = Upl{N, A)U for N even. (12) 

Ps{N,-A) = (a^)®"J7p:(7V,A)t/(a^)®" for odd 

(13) 

where U = cr^ / ® cr^ (g) .... is a tensor product over all 
odd sites, which has the property UH{A)U = —H{—A), 
and p* a complex conjugated matrix in the basis where 
is diagonal. Note that the above transformations depend on 
parity of N and, in particular, result in the current being an 
odd (even) function of A for odd (even) N, 



j(A) = (-l)^j(-A) 



(14) 



For the case = 3 the system of ^3 = 32 algebraic equa- 
tions obtained with the help of the symmetry ^ is reported 
in Appendix A. From the solution of this system one can eas- 
ily check that the analytical expression for the average current 



J(7, A) = Tr{p{a^{i)ay{i + 1) - <Ty{i)a^{i + 1))) is of the 
form 



J(7,A) 



1=0 ' 



n(A)7 



2(i+l) 



A^ = 3, 



(15) 



with coefficients ai{A), f3i{A), odd and even polynomials in 
A, respectively. This implies that the current is an odd func- 
tion of A, in accordance to (fT4l) . The first leading coefficients 
relevant for strong and weak coupling limits are 

ao(A) = 144A5(A2 - 1)\9A^ + 59), 

ai(A) = 8A^{A^ ~ 1)2(9420 + 13895A2 + 5822A'^ 

+6561A^ + 2396A^ + ISGA^^), 



a9(A) = 248832A(265 + 96A2), 
aio(A) ^ 8957952A. 



(16) 



and 



;3o(A) = 2AA\A^ - 1)'*(416 + 155A^ + 9A^), 
Pi{A) = 2A^{A^ - 1)2(34032 + 45152A2 - 7764A'^ 



Ao(A) = 62208(228A4 + 644A2 + 355), 
Ai(A) = 2239488(1 + 2A2). 



(17) 



To avoid lengthy expressions we omitted other coefficients 
since they can be easily obtained from the exact solution of 
the algebraic system in Appendix A. Substituting the leading 
coefficients into Eq. ( fTSl l we obtain expressions for the current 



J(A,7) 
J(A,7) 



472A[18(9A2 + 59)(A - A^)^ + j^Fi] 



12A2(A2 - 1)2(9A4 + 155A2 + 416) 

4A(96A2 + 3672 + 265) 
228A4 + 3672(1 + 2A2) + 664A2 + 355 



72^2' 
(18) 



valid for arbitrary values of A and small, 7^1, and large, 
7^1, values of 7, respectively. Here Fi = 186A^° + 
2396A*^ + 6561A6 + 5822A4 + 13895A2 + 9420 and = 
153A12 + 3411A10 + 19179A« + 24007A6 + 7764A'* + 
45152A2 + 34032. Note, from the second Eq. in that in 
the infinite coupling limit 7 00 the current is: 



lim J(7,A) 

7—^00 



4A 



1 + 2A2 ^ 



(19) 



this coinciding with the expression derived in 

IM by means 

of a perturbative approach. Similar results are obtained for 
other values of A^, with the only difference that the case of A^ 
even the coefficients of the series in Eq. ( fTSl l are even func- 
tions of A, e.g. for even (odd) A^ the current is an even (odd) 
function of A. Explicit expressions of p (and derived quan- 
tities), however, become very involved as A^ increases since 
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Figure 3: (Color online) Top panels. Stationary spin currents J (left 
panel) and von Neumann entropies (right panel) vs A for the "hLxxz 
chain of length A*' = 4 for different values of 7 indicated in the fig- 
ure. Curves are obtained from exact stationary solution of the Lind- 
blad Master Equation. Bottom panels. Corresponding von Neumann 
entropies vs stationary currents J for different coupling constants 
7 indicated in the panels. Bottom right panel show details of the 
crossover through the cusp point at which minimum of VNE and 
maximum of current are in exact coincidence (e.g. they occur at the 
same value of A). All curves are obtained from exact stationary so- 
lutions of the Lindblad equation. 



the number of algebraic equations to solve increases expo- 
nentially (see Eq. (|9]l. Although it may be a problem to get 
analytical expressions for such large systems, it is not a prob- 
lem to solve them numerically with very high accuracy to be 
considered exact for any practical purpose. 



III. TRANSPORT PROPERTIES AND VON NEUMANN 
ENTROPY FOR OUT OF EQUILIBRIUM FINITE SIZE XXZ 
CHAINS 

The symmetry approach discussed above allows to obtain 
exact analytical results for the physical quantities of interest 
in the transport, such as the spin current J and the von Neu- 
mann entropy VNE{pst) = — l^i-^j' where Ai denote 
the eigenvalues of the stationary reduced density matrix. 



In the top panels of Fig. [T] we show stationary spin cur- 
rent J (left) and von Neumann entropy (right), as obtained 
from the exact solution of the system in Eq. ( IA3l l for the case 




A 



Figure 4: (Color online) Top panels. Stationary spin currents J 
(left panel) and von Neumann entropies (right panel) vs A for the 
T-Lxxz chain of length N — 5 for different values of 7 as indi- 
cated in the figure. Curves are obtained from exact stationary so- 
lutions of the Lindblad Master Equation and, taking into account 
their antisymmetry property in A, they have been plotted for pos- 
itive A values only, just for graphical convenience). Bottom Pan- 
els. Details of the reversal current transition (left panel) and cor- 
responding von Neumann entropy. Curves from left to right in the 
left panel refer to 7 values 1.9, 1.95, 1.96285, 1.965, 1.97, respec- 
tively. Curves from top to bottom in the right panel refer to 7 values 
1.8, 1.9, 1.97, 2, 2.5, 3, 5, 6, respectively. 



= 3, as function functions of the parameters 7, A. Sections 
of the spin current J (left) and von Neumann entropy (right) 
surfaces are reported in the bottom panels for different values 
of parameter 7. We see that the mean current is an odd func- 
tion of A and the VNE curve has minima very close (practi- 
cally in correspondence) to maxima of | J|. This is clear also 
from Fig. |2] where a parametric curve J vs VNE is reported 
for different values of 7. Notice that in the strong coupling 
limit the VNE becomes zero exactly in correspondence to the 
maxima of | J|. This implies that optimal transport (e.g. the 
maximal current) is achieved in correspondence of pure states 
of the open XXZ chain . Notice that the current vanishes for 
A = ±00 (Ising limit) and the VNE attains its maximum at 
A = (free fermion point). It is also worth to note here that 
although analytic results are available for the Lindblad Master 
equation at the free fermion point lfl9l - l23ll , these results do not 
apply to our choice of Lindblad operators (O. 
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Figure 5: (Color online) The von Neumann entropy VNE vs sta- 
tionary current J for the open T-Lxxz chain of length N — 5 with 
different coupling constants 7 as indicated in the panels. The curves 
are obtained as parametric plots for fixed 7 and A varied in the range 
] — cxD, cxj [. From panel insets one can see that the minima of the VNE 
are always very close to maxima of the absolute value of the current. 
As the coupling constant 7 is increased the values of the VNE at 
the minima decrease and in the strong coupling limit 7 — >■ cxj they 
reduce exactly to zero (pure states). 



Similar results, obtained for the case = 4, are shown 
in Fig. [3] We see that the spin currents J (top left panel) 
for different values of 7 are even functions of A, and that the 
von Neumann entropies (top right panel) attain their maximal 
value at A = ±00 where J is zero, while their minima are 
still very close to the maxima of | J|. Notice from the bottom 
left panel, that while the minima of the VNE always reduce to 
zero in the strong coupling limit 7 00, the correspondence 
between the max of | J| and the min of VNE remains non ex- 
act also in this limit. This property is completely general for 
all > 3 (see below). It is also worth to note from this Fig. 
that as 7 is increased from small (green curve) to strong cou- 
pling Umit (red curve) the VNE- J curve undergoes a crossing 
point with the development of a cusp at a particular value of 
A for which the maximal current exactly coincides with the 
minimum of the VNE. This process is clearly shown in the 
bottom right panel of Fig. [3] where the cusp appears on the 
black curve 7 = 0.85 for J w 0.6. This folding of the curve 
as 7 is increased is also a general property observed for other 
values of N. 



Figure 6: (Color online) Top panels.Stationary spin currents J (left 
panel) and von Neumann entropies (right panel) vs A for the Hxxz 
chain of length AT = 6 for different values of 7. Curves are obtained 
from exact stationary solution of the Lindblad Master Equation. Due 
to the parity of J( A) for A'' even, we have restricted the figure to the 
range A > only. 



As A^ is increased, however, curves becomes more and 
more complicated with the occurrence of the interesting phe- 
nomenon of current reversal. This is first observed for the case 
N — 5 reported in Figs. |4]|5]where the spin current J and the 
von-Neumann entropy VNE, are depicted both as explicit and 
as implicit functions of A. From the top left panel of Fig. 2] 
we see that, quite interestingly, the system undergoes current 
reversals in the critical region — 1 < A < 1 (see the reversal 
occurring at A w 0.4), this increasing by two the number of 
peaks of | J| discussed for the A^ = 3, 4 cases (notice that, due 
to the antisymmetry of J(A) for A^ odds, we restricted these 
figures to the range A > 0). The negative peaks in the current 
correspond to a relative minima appearing in the VNE curve, 
as one can see from the top right panel of this figure. Notice, 
however, that the correspondence between secondary peaks 
of |J| and corresponding VNE minima occurs with delay in 7 
(notice that there exist tiny minima in J which have no cor- 
responding minima in the VNE curve). Details of the current 
reversal phase transition and the behavior of the correspond- 
ing VNE curves are reported in the bottom panels of Fig. |4]for 
different values of the coupling parameters. Also, in Fig. |5]we 
have depicted the VNE vs stationary current J as parametric 
plots for fixed 7 and A varied in the range ] — 00, oo[. From 
the panel insets one can see that the minima of the VNE are 
always very close to maxima of the absolute value of the cur- 
rent. As the coupling constant 7 is increased the values of the 
VNE at the minima decrease and in the strong coupling limit 
7 — > 00 they reduce exactly to zero (pure states), in analogy 
with the cases A^ = 3, 4 discussed before. 

Similar behavior, with appearance of reversal currents in 
the critical region is observed also for the case A^ = 6, as one 
can see from Fig. |6] The appearance of the reversal current 
can be seen as a precursor of a quantum phase transition driven 
by the boundaries as the coupling parameter is varied. As a 
general rule we conjecture that the number of extrema of | J| 
is at most N ~ 2 times for A^ even and A^ — 1 times for N 
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netic reservoirs by means of the Lindblad master equation. 
Exact solutions for the stationary density matrix have been 
constructed for chains of small sizes using the quantum sym- 
metry of the system. Using these solutions we have investi- 
gated the transport property of the chain in terms of the von 
Neumann entropy. As a result we have shown that the maxi- 
mal spin current always occurs in the proximity of a minimum 
of the VNE and for particular choices of the anisotropy param- 
eters it can exactly coincide with it. More precisely, we found 
that as the coupling constant increases, the VNE- J curve un- 
dergoes a crossing point with the development of a cusp for a 
particular value of A for which the maximal current is in exact 
coincidence with the minimum of the VNE. We also showed 
the existence of current reversals e.g. the presence of nega- 
tive peaks in current, which correlate to relative minima in the 
VNE curve, when the anisotropy parameter is in the critical 
region —1 < A < 1. These relative minima disappear in 
the small coupling limit, while in the infinite coupling limit 
we show that the minima of the VNE becomes exactly zero, 
meaning that maximal transport in this case is achieved with 
states very close to pure states. 



Figure 7: (Color online) von Neumann entropy VNE vs stationary 
current J for the open T-Lxxz chain of length A'^ = 6 for coupling 
constants 7 as indicated in the panels. As the coupling constant 7 is 
increased the values of the VNE at the two minima decrease and in 
the strong coupling limit 7 — >■ 00 they reduce exactly to zero (pure 
states). 



odd. These extrema occur always in the critical region and 
may disappear as the coupling constant is decreased. 
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Appendix A: 

Exact stationary density matrix elements for the case A'^ = 3 

One can check that the most general form of the density matrix compatible with the symmetry operation T in Eq. (|4]i for the 
case N = 3,is 
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aii 


ai2 + ibi2 


ai3 + ibiz 


014 + ibiA 


015 - 


1- ibi5 


016 + ibi6 


Ol7 


+ ibn 


(1- 


i)ai8 \ 


ai2 


- ibi2 


a22 


023 + *^23 


(1 + i)024 


025 " 


t- ib25 


026 + ib26 


027 


+ ib27 


ibiA 


- 014 


ai3 


- ibis 


0-23 - ib23 


033 


1026 + ^26 


035 - 


'r ib35 


(1 + 2)036 


037 


+ ib37 


ibiQ 


- 016 


ai4 


- ibii 


(1 - z)a24 


626 - ia2(, 


022 


627 - 


- «027 


023 + ^^23 


025 


+ ib2b 


mi2 


+ 612 


ai5 


- ibiz 


025 - ib2b 


035 - «^35 


1027 + ^27 


055 


«037 + 637 


(1- 


h i)057 


ibn 


- Oi7 
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- ibie 


026 - ib2Q 


(1 - z)a36 


023 - ib23 
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- 1037 


O33 


035 


+ ib35 


iai3 


+ bi3 


ai7 


- ibn 


027 — ib27 


037 - 1^37 


025 - ib2b 


(1- 
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035 - ib35 




055 


iai5 


+ bi5 


V(i- 


h i)ai8 


—an — ibu 


-ai6 - «foi6 


bi2 - iai2 


-Ol7 


- i6i7 


bi3 ~ iai3 


bi5 


- »Oi5 


On / 



(Al) 



This can be easily proved by substituting an arbitrary form of p, e.g. with arbitrary diagonal real elements, a^, and arbitrary off 
diagonal complex elements, Oy + ibij, into Eq. ^ and by imposing the equality. This gives the set of relations reported in Eqs. 
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([Tol l. (fTTT i and made explicit in the above p. 

The exact expressions of pij are obtained by substituting (Al) into the Lindblad master equation and by imposing the station- 
arity. This gives the following system of algebraic equations : 

2(ai4 - an) + 2{bu + 617) + 027 - 627 + 6ai8 = 0, 
2aii - a22 - a55 = 0, 
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(2a22 - 


- On - 


- 033) + 4^23 = 0, 
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,2055 - 


- On - 


- 033) - 453,5 = 0, 
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[2a33 - 


- 022 - 


-055) + 4(635 -^23) = 


= 0, 
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2a24 - 


- Oi3 - 


- ^13) + 2(626 - 026) = 


= 0, 
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,2057 - 


- Oi3 - 


- ^13) + 2(037 - ^37) = 


= 0, 








7 


(ai2 + 


615 - 


5)16 - 3625 + ^26 + ^37 


) + 023 - 


- 035 = 


= 0, 




7 


(ai5 + 


026 - 


037 - 3oi6 - 612) - 614 - bn 


+ 2A6 


16 = 0, 




7 


(3ai5 - 


- 026) 


+ 2(6i3- A615) = 0, 










7 


(36l2 - 


- 037) - 


- 2(ai3 - Aai2) 0, 










7 


(3a37 - 


- bi2) - 


- 2(057 - 036 + 627 - 


A637) = 


= 0, 






7 


(3a26 - 


- O15) 


+ 2(024 - 036 + ^27 - 


A526) = 


= 0, 






7 


(3a25 - 


- O15 - 


- 026 — O37 — 612) — 635 + 623 


= 0, 






7 


(2ai3 - 


- O24 - 


- 057) + 2(5i2 + 5i5 - 


2A613) 


= 0, 






7 


(26i3 - 


- 024 - 


O57) - 2(oi2 + Oi5 - 


2A013) 


= 0, 






7 


2(a24 


+ 013) 


- 4oi4 - Oi7 - O23] - 


- 2(6i6 - 


- Abu) = 0, 




7 


2(a24 


+ 613) 


- 45i4 + bn - b2z] + 


2(oi6 - 


A014) 


= 0, 




7 


2(057 


+ 013) 


- 46i7 + bu - 635] + 


2(oi6 - 


A017) 


= 0, 




7 


2(057 


+ 613) 


+ 4oi7 + Ol4 - 035] + 


2(6i6 - 


A617) 


= 0, 




7 


2(057 


+ 613) 


- 4o35 + 023 + 017] 4 


2(625 - 


A635) 


-0, 




7 


2(057 


+ 013) 


- 4^35 + 623 - bn] + 


2(033 - 


055 - 


025 + A035) 


7 


2(024 


+ 013) 


- 4023 - Oi4 + O35] - 


- 2(625 - 


- A623 


) = o, 




7 


2(024 


+ &I3) 


- 4523 - bli + 635] + 


2(022 - 


033 + 


025 - A023) 


7 


2(011 


+ 055) 


- 55i5 + 626] + 2(oi3 


- A015 


)-o, 






7 


2(oii 


+ 022) 


- 5oi2 + 637] - 2(613 


- A5i2 


)-o, 






7 


2(033 


+ 022) 


- 5526 + 5i5] + 2(024 


+ 027 - 


036 - 


A026) = 


0, 


7 


2(033 


+ 055) 


- 5637 + O12] + 2(036 


- 027 - 


' 057 4- 


A037) = 


0, 


7 


2(oi4 


- 017) 


+ 6027 + O18 - 036 + 


2(623 + 


^'35)] H 


h 2(626 - 


^37 


7 


2(023 


+ 035) 


- 6627 + O18 - 036 + 


2(6i4 + 


bn)] 4 


- 2(026 - 


037 


7 


2(023 


+ 035) 


- 6036 + 027 - b27 + 


2(623 + 


635)] - 


- 2(026 " 


037 


7 


(ai2 + 


615 - 


36l6 - ^25 + &26 + ^37 


) + 014 - 


f 017 - 


- 2A016 = 


= 0, 



(A2) 



0, 



for the real and imaginary parts, aij,bij, respectively, of 
the density matrix elements. One can easily check that the 
solution becomes unique when the normalization condition 
Tr{p) = 1 is imposed, this completely solving the problem. 
This approach, being general, can be used for other values of 
N. 



1611). in the form 



P7- 



i/+i(/-a^)^ 



1 + 2A- 



:(a"-a^)(»(/+a"). (3) 



From this expression we readily compute the eigenvalues 
for the reduced density matrix, = (0, 0, 0, 0, 0, 0, Ai, A2) 
where 



In the limit 7 — )• 00 the above set of equation simplifies 
and yields a compact solution for the density matrix, (see also 



Al,: 



1 + 4A2 + 4A4 ± 2A^/2Vl + 4A2 + 4A4 
2(1 + 4A2 +4A4) ■ 



(4) 



8 



At the point A = l/\/2, Ai = 0,A2 = 1 and the density 
matrix spectrum becomes that of a pure state. The point A = 
1 / \pl is exactly the point where the steady current attains its 
maximum value 7Tiaa;o<A<oo^( A) = \f2. 

Analogous calculations for higher number of sites > 3 
lead to similar results, with the difference that the steady 
state in the limit of large couplings becomes pure not ex- 
actly at the point A where the magnetization current has an 
extremum, but quite close to it. For A^ = 4, the extremum 
of the current at 7 oo reaches its maximum at the point 



A = (-1/4 + 73/2)1/2 ^ Q ^jjjjg jjjg density matrix 
spectrum becomes "pure" at A* = -\/3/2 ~ 0.866. Note how- 
ever that the maximum peak J( A) is rather broad in this case, 
so that the maximal j value and its value at the "pure" point 
J(A*) are very close, (Jmax- J(A*))/Jmax « 0.009. At 
A^ = 5, there are two extrema, at Ai_2 ~ 0.384, 0.814, while 
the state becomes pure at the points A* 2 ~ 0.383, 0.924. In 
the point Ai w 0.384 the peak is narrow, and so Ai w A*, 
while the second peak at A2 ~ 0.814 is broad one, justifying 
the discrepancy between the Ai and Aj (data not shown). 
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